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We investigate non-equilibrium nature ol fluctuations ol black hole horizons by applying the 
fluctuation theorems and the Jarzynski equality developed in the non-equilibrium statistical physics. 
These theorems applied to space-times with black hole horizons lead to the generalized second law 
of thermodynamics. It is also suggested that the second law should be violated microscopically so 
as to satisfy the Jarzynski equality. 
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Introduction. - There arc striking similarities be- 
tween space-time with horizons and thermodynamic sys- 
tems. The analogies have been extensively investigated, 
especially, in the black hole thermodynamics. A black 
hole behaves as a blackbody with the Hawking tem- 
perature Th = fin /2tt [l|, and energy AE flowing into 
the black hole can be identified as the entropy increase 
ThASbh of the black hole. Here, K is the surface grav- 
ity at the horizon and the entropy of the black hole 
Sbh is proportional to the area of the horizon A as 
Sbh = A/AGh in the Einstein-Hilbcrt theory of gravity. 
The thermal behavior is essentially quantum mechanical. 
Furthermore, such a thermal behavior is not restricted to 
globally-defined horizons like an event horizon of a black 
hole, but also applicable to local horizons such as the 
Rindler horizon of a uniformly accelerated observer. The 
point is emphasized in e.g. Q. At the classical level, 
this is related to an interplay of a local change of the 
horizon area and the Einstein equation. In particular, 
T. Jacobson pointed out Q that the Einstein equation 
is derivable from a thermodynamic relation for the local 
Rindler horizon. At the quantum level, the local notion 
of horizon entropy has more fundamental meaning, since 
it gives transition rates of area-changing irreversible pro- 
cesses of black holes, and will be related to the quantum 
statistical nature of the space-time. If there exists a fun- 
damental quantum structure of space-time behind the 
horizon thermodynamics, how can we probe such micro- 
scopic states? 

In order to tackle this question, let us go back into history 
at the early twenties century. Einstein proposed several 
methods to determine the Avogadro constant Na- In 
particular he investigated the dynamics of the Brown- 
ian motion in His aim was giving a proof of the 
reality of atoms. If we knew nothing about the atoms, 
how could we count the number of atoms contained in 
one mole gas? His idea was to measure fluctuations of 
physical quantities around thermal equilibrium. In ther- 
modynamics, the unit of energy per kelvin per mole is 



given by the gas constant R, and any thermodynamic 
observable in equilibrium is expressed in terms of R. As 
far as one measures thermal averages in equilibrium, it is 
impossible to disentangle the macroscopic unit of energy 
R into the Avogadro constant Na and the Boltzmann 
constant ks- Einstein's insight was in paying attention 
to fluctuations around thermal equilibrium. A thermal 
quantity fluctuates because of the microscopic structures 
behind the thermodynamics. He reversed the entropy 
formula of Boltzmann S = ks log W as W — exp(S , /fcs), 
and obtained the mean squared values of fluctuations in 
terms of the microscopic unit of energy ks and various 
response functions For example, the entropy fluctu- 
ation at constant pressure is given by ((AS) 2 ) = krsCp 
where C p is the heat capacity at constant pressure. This 
is now understood as a fluctuation-dissipation theorem, 
which relates the fluctuation around the equilibrium to a 
response function. 

After Einstein, many important works have been per- 
formed to understand out-of-equilibrium properties of 
fluctuations. One of the most notable developments will 
be the non-equilibrium fluctuation theorem first discov- 
ered in numerical simulations @ . There are several vari- 
ations of the theorem. An adequate one in the present 
context is the Crooks fluctuation theorem Q- Assume 
that the system is initially in thermal equilibrium at in- 
verse temperature (3 with an external parameter X F (0). 
Then one changes the external parameter A F (t) as a func- 
tion of time from t = to t = T. The system evolves 
according to the Hamiltonian dynamics. The procedure 
of changing the parameter corresponds, for example, to a 
process of moving a piston and it needs not to be quasi- 
static. In changing the external parameter, the system 
becomes out-of-equilibrium. For each microscopic state, 
one measures the amount of work exerted on the system, 
and takes an ensemble average over the initial density 
matrix. Define p F (W) as a probability that the exerted 
work is given by W under the change of parameter A F (t) . 
X F (t) is called a forward protocol. We also consider a re- 
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versed protocol in which we change the external parame- 
ter in a reversed way as A R (i) = A F (T — t) from t = to 
t = T. The system is assumed to be initially in thermal 
equilibrium at the same temperature, but with a different 
external parameter A R (0) = A F (T). Define /5 R (— W) as a 
probability that the work is given by — W in the reversed 
protocol. The Crooks fluctuation theorem states that the 
ratio of these two probabilities is given in terms of the 
work and the difference of free energies F(A) between the 
two equilibrium states, 



P F (W) 
p R (-W) 



J(W-AF) 



(1) 



where we defined AF = F(A R (0)) - F(A F (0)). From 
the Crooks fluctuation theorem, the Jarzynski equality 
H can be obtained, (e~P w ) = e~@ AF . Here, the an- 
gled bracket stands for the average with the probability 
p F (W). It is surprising since the average of exponenti- 
ated work in non-equilibrium processes in the left hand 
side is related to the difference of equilibrium quantities 
at the beginnings of the protocols. These relations imply 
that non-equilibrium fluctuations are miraculously ar- 
ranged so as to satisfy such very nontrivial relations. By 
using the Jensen's inequality, (cxp(x)) > exp((x)), the 
Jarzynski relation is reduced to ((W — AF)) > 0, which 
implies the second law of thermodynamics. Though the 
average of W — AF is always non- negative, the Jarzynski 
relation requires that there is a nonzero probability for 
the quantity to take a negative value microscopically. 

The purpose of this letter is to apply these celebrated 
relations of non-equilibrium fluctuations to space-times 
with horizons. In particular, we consider non-equilibrium 
fluctuations of black hole horizons. 

Transition Rates. — We consider a coupled system of 
a black hole and matter. The external parameter A(t) 
characterizing the system Hamiltonian, which appeared 
in the fluctuation theorem, can be arbitrarily chosen 
here such as a height of a potential for the matter. If 
the whole system is controlled by a unitary time evolu- 
tion with time-reversal symmetry, a transition probabil- 
ity VVxF(t) (C — > C) from one configuration C to another C 
under a time evolution of the external parameter A F (t) is 
equal to a probability W^rm (C — > C) from C to C under 
the reversed change of the parameter X K (t). In the pres- 
ence of a black hole, however, the time-reversal symmetry 
is violated by imposing the ingoing boundary condition 
at the horizon. In a black hole space-time, regular coor- 
dinates near the horizon, i.e. Kruskal coordinates (U, V), 
are defined by U = -K,- 1 e-^ t ~ r ^ and V = n^e^+^l 
Here, t and are the Schwarzschild-time and the tortoise 
coordinates. Quantum fields near the horizon are classi- 
fied into two types of chiral fields (in a two-dimensional 
sense on (t, r*) plane), one depending on U and the other 
on V. Fields depending on V are ingoing waves falling 



into the black hole while those depending on U are propa- 
gating nearly along the horizon and correspond to outgo- 
ing modes. The regularity at the horizon requires occu- 



pation of outgoing modes 4>(u>) 



-iojU 



to vanish at the 



(future) horizon. Namely, we must impose the vacuum 
condition for the outgoing modes in the Kruskal coordi- 
nates. On the contrary, there is no constraint for the in- 
going modes, and the conditions are asymmetric between 
U and V. The time-reversal transformation t —¥ —t ex- 
changes the coordinates U and V, and the presence of 
horizon violates the time-reversal symmetry of the quan- 
tum states. Therefore, the above transition probabilities 
are not necessary the same. 

The ratio of the above transition probabilities was evalu- 
ated by Massar and Parentani [§] for Hawking radiation 
processes. They have shown that the transition rates 
for systems with a black hole horizon are governed by 
changes in the horizon area. In the present case, the ra- 
tio of a transition probability between a configuration C 
with black hole area A and another one C with area A' 
under a fixed value of external parameter A is given by 

W X (C(A) -> C'(A')) 



= cxp(AA/4Gh), 



(2) 



Wx{C'{A')^C{A)) 

where the change of area AA = A' — A is assumed to 
be small. In deriving this, they used the WKB approxi- 
mation for the system wave function, and calculated the 
transition rates in the first Born approximation for the 
interaction between the detector and radiation field. A 
similar result was obtained in a different way in (lfjj . If we 
identify AA as the energy AF emanated from the black 
hole by AA/AGh = —AE/T H , it becomes the Boltzmann 
factor exp(—AE/Tjj) in the Hawking radiation. The ra- 
tio eq.© takes into account the back reaction of the ge- 
ometry against emanating radiation quanta. If detailed 
balance is satisfied in the processes, the ratio eq.® is 
identified as the ratio of a probability in the configura- 
tion C'(A') to that in C(A), and hence consistent with 
the entropy of a black hole Sbh = A/4GH. 

In the proof Q , they have used an observation by Carlip 
and Teitelboim [ll| that, if we consider a coupled system 
of a black hole and matter exchanging energy between 
them, one needs to add a boundary term to the bulk ac- 
tion, S = Sbuik + AQ/8ttG. Here, 6 = nt for on-shell and 
stationary metrics. Then in quantizing the system, the 
Wheeler-DeWitt equation H to t^ = in the bulk must 
be supplemented by the boundary Schroedinger equation 



111 ] ihd&'fy = — (A/8irG)^ , and the total system's wave 
function evolves as exp(iAQ/8TrGh) in the WKB approx- 
imation. 

Here, we argue that the ratio ([2]) should be valid also for 
processes including classical absorption of energy into the 
black hole, if we generalize the notion of transition prob- 
abilities in the following way. The matter system outside 
the horizon dissipates the energy by transferring it into 
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the black hole. Furthermore, the system feels a kind of 
thermal noise due to the Hawking radiation. Hence, by 
including both effects of the heat transfer, in and out, at 
the horizon, the effective equation of motion for matter 
is controlled by a stochastic Langevin equation with dis- 
sipation and noise terms. In such a situation, one can 
define a probability distribution of the system to take 
some configuration. Time evolution of the probability 
distribution function is described by the Fokker-Planck 
equation. Clearly the time reversal symmetry is violated, 
and we can expect an asymmetry between the probabili- 
ties of the forward and the reversed processes. The ratio 
is evaluated in general Langevin processes in [12j. By 



applying it to our case, the energy transfer into the black 
hole can be rewritten as the area change of the black 
hole through the first law of black hole thermodynam- 
ics (ASbh = AE/Tbh)- It can be possible to bring 
it together with Hawking radiation effect into the total 
amount of area change. This suggests the validity of the 
probability ratio ([2]) for wider situations including classi- 
cal absorption of energy into the black hole. Details will 
be discussed in a forthcoming paper. 

N on- equilibrium Fluctuations of Horizons. — We con- 
sider a sequence of configurations of a coupled system 
of a black hole and matter, and denote it as T = 
{C (A ),Ci(Ai), . . .,C M (A M )}- The configuration C k is 
realized at a discrctized time t — kAt. Each transition 
probability is given by W A p( tje ) {C k {A k ) ->• C k +i(A k+1 )). 
If we assume the Markov process, the transition proba- 
bility for the sequence of configurations T to be realized 
is given by a product of them, 

M-1 

P F (D = J] W XF(tk} (C k (A k ) C k+1 (A k+1 )). (3) 

k=0 

As we argued, the sequence can represent a general 
process of absorbing and emitting matter through the 
black hole horizon. On the other hand, the probabil- 
ity for the reversed sequence of configurations T* = 
{Cm{Am), ■ ■ ■ ,Ci(Ai),C (Ao)} with the reserved change 
of the external parameter is given by 



A/-1 



P R (r*)= [] W x n {tk] (C M _ k (A M _ k ) ^ C M _ k+1 (A M - k+1 )) 

k=0 
M-1 

[ e^W xF{tk) (C k+1 (Ak+i) ~^ C k (A k ))). 



k=0 



(4) 



Here we have used eq. @- Hence the ratio of these two 
probabilities is given by 

p p (r) ( a m -Aq \ _ 

= CX P I — ZFTE — ) = exp(5p(r)). (5) 



p R (r*) 



4Gh 



Sp(T) is defined as the logarithm of the ratio, and pro- 
portional to the difference of area, which is not necessar- 
ily small. 

We now derive a Crook's type fluctuation theorem in the 
black hole system. The matter is assumed to be in ther- 
mal equilibrium with Hawking temperature Th with an 
external parameter A^O) at the beginning. First define 
the total dissipation AS(T) by 

«p(-AS(T)) = exp(-S P (T)), (6) 

where p_x(p, R-)(t ) is the initial probability distribution for 
matters under the forward or reversed protocols. We 
assume that these probability distributions are canonical 
distributions with the Hawking temperature. A relation 
A5(r*) = -AS(T) is followed by S P (T*) = -S P (T). 
A transition probability to produce the total dissipation 
AS(T) = AS under the forward protocol A F (i) is now 
given by 

P F (AS) = 5> AF(0) (C )P F (r)«5(AS(r) - AS)). (7) 
c ,r 

Here, X)c P\ F (o)(fio) stands for a sum over all possible 
initial states weighted by the initial distribution. is 
a path integral for all possible trajectories. By using (O 
and ([5]), it is straightforward to show the Crook's type 
fluctuation theorem 



p F (AS)e 



-AS _ „R 



P«(-AS). 



(8) 



From the definition of AS(T) in ([6j), it is identified as a 
sum of the black hole entropy and the entropy of matter; 



AS(T) 



AA 

4GH 



(3 H {AE - AF). 



(9) 



By integrating the equation © over AS, it gives a 
Jarzynski type equality 



<e-**> = 1. 



(10) 



The Jensen inequality implies the generalized second 
law of thermodynamics (AS) > of the black hole 
with matter [l3|. An important point here is that it is 
satisfied only in an averaged sense, and in order to satisfy 
the Jarzynski type equality (|10[) . entropy decreasing 
processes must exist as individual processes, and their 
probabilities are miraculously arranged to satisfy the 
Jarzynski equality. 

Conclusions and Discussions. - In this letter, we 
have applied the recent developments in non-equilibrium 
statistical physics to area changing processes of a 
black hole interacting with external matter. The 
non-equilibrium fluctuation theorems a la Crooks and 
Jarzynski are derived, which lead to the generalized 
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second law of black holes. The second law holds only 
after taking a thermodynamic average, and it should be 
violated as an individual process in a miraculous way to 
satisfy the Jarzynski equality. 

There are several issues that further investigations are 
necessary. First, we have applied the formula ([2j to 
general processes including absorption and emission 
of classical and quantum particles. As argued, it will 
be proved by using a path integral formulation of 
stochastic processes for the matter fields in a black 
hole background. More details will be reported in a 
separate paper. Instead of using the stochastic equation 
for the matter that may be obtained after integrating 
out internal degrees of freedom of a black hole, we may 
consider a whole quantum system containing both of the 
internal and external degrees of freedom. This will be 
achieved by considering the global Hilbcrt space of the 
left and right wedges altogether. 

Another issue is to apply a different type of fluctuation 
theorems, such as the steady state fluctuation theorem. 
If the matter continues to be in a thermal equilibrium 
with a different temperature from Th, there is a constant 
flow of energy between the black hole and the matter. 
The fluctuation theorem is also applicable to such a 
situation. From the steady state fluctuation theorem, we 
can derive a fluctuation-dissipation relation by assuming 
Gaussian fluctuations around a long-time average in a 
late time. We may also be able to obtain a Grecn-Kubo 
relation that relates the correlation of energy flow and 
a proportionality coefficient between area change of 
horizon and energy flow. 

More ambitiously, as Einstein asked himself how to 
prove the reality of atoms, we may also ask ourselves, 
what is an analog of the Avogadro constant in the 
gravity, or how we can probe such quantities by looking 
at some kinds of fluctuations caused by the microscopic 
structures of the space-time. An analogous quantity to 
the Avogadro constant will be the black hole entropy per 
unit area ~ 10 70 [l/m 2 ] (or its exponential). Since it is a 
huge number, the fluctuation will be largely suppressed. 
Finally, note that the fluctuation of the black hole 



entropy (horizon area) is proportional to the specific 
heat, which is negative for the Schwarzschild black hole. 
As it means an instability of the black hole, it becomes 
the more important to take into account non-equilibrium 
effects of the fluctuations. 
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